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For the spin-glass chain in an external field h, a non-zero weight at the origin of the bond 
distribution p(J) is known to induce a non-analytical magnetization at zero temperature : for 
p(J) ~ yl|J|^~^ near J ^ 0, the magnetization follows the Chen-Ma scaling M ~ /i'^/(2+m) ^j^jg 
paper, we numerically revisit this model to obtain detailed statistical information on the ground 
state configuration and on the low-energy two-level excitations that govern the low temperature 
properties. The ground state consists of long unfrustrated intervals separated by weak frustrated 
bonds: We accordingly compute the strength distribution of these frustrated bonds, as well as the 
length- and magnetization- distributions of the unfrustrated intervals. We find that the low-energy 
excitations are of two types (i) one frustrated bond of the ground state may have two positions 
that are nearly degenerate in energy (ii) two neighboring frustrated bonds of the ground state may 
be annihilated or created with nearly zero energy cost. For each excitation type, we compute its 
probability density as a function of its length. We moreover show that the contributions of these 
excitations to various observables (specific heat, Edwards- Anderson order parameter, susceptibility) 
are in full agreement with direct transfer matrix evaluations at low temperature. Finally, following 
the special bimodal case ±J, where a Ma-Dasgupta RG procedure has been previously used to 
compute explicitly the above observables, we discuss the possibility of an extended RG procedure : 
we find that the ground state can be seen as the result of a hierarchical 'fragmentation' procedure 
that we describe. 

I. INTRODUCTION 

In this paper, we consider the one dimensional spin-glass chain in a small external field h > 

H = -'^JiaiUi+i- h'^Gi (1) 

i i 

to obtain detailed results on the ground state and the low-energy excitations, as a function of the exponent /i > 
characterizing the weight of the coupling distribution for small couplings 

(2) 

As is well-known, the previous model is equivalent to a random-bond and random-field ferromagnetic chain 0| 

H = -'^\Ji\SiSi+i - h'^XiSi (3) 

i i 

where Xi = 11^=1 sgii( J,). 

A. Bimodal distribution Ji = ±J: Imry-Ma argument and real-space RG 

For the special case of the bimodal distribution Ji — ±J with probabilities (1/2, 1/2), the model Q corresponds 
to a pure Ising chain | J^j = J in a bimodal random field hi = hxi = ±h. The Imry-Ma argument Q for the random 
field Ising chain can be immediately translated for the spin-glass in external field, since the domain walls of the RFIM 
now becomes frustrated bonds for the spin-glass : the random magnetization m of an unfrustrated domain of length 
I is of order m ~ Vl, i.e. it gives rise to an energy of order 2h\/l in the external field h, whereas a pair of two 
frustrated bonds has for energy cost 4 J. As a consequence, the ground state is made of unfrustrated domains having 
the typical Imry-Ma length Ljm ~ 4 J'^ /h?. The real-space Ma-Dasgupta RG 2J allows to construct explicitly the 
positions of frustrated bonds and to compute various statistical properties, such as the distribution of the domain 
lengths. This approach moreover yields the statistics of low-energy two- level excitations Q. A natural question is 
then : what are the corresponding results for a general distribution p{J) that is not bimodal ? It turns out that a 
different behavior occurs if p{J) has some weight at small couplings J ~ 0. This case, which includes the Gaussian 
distribution, completely changes the physics of the model, as we now discuss. 
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B. Distributions with small couplings p(J) — yl| J]'' : Chen-Ma argument 

For distributions presenting some weight at small couplings (01 the above Imry-Ma argument for the bimodal case 
is replaced by the following Chen-Ma argument Q. The essential idea is that frustrated bonds will be now located 
on weak bonds, in contrast with the bimodal case where the cost of a frustrated bond is the same everywhere. More 
precisely, the Chen-Ma (CM) argument is as follows : the bonds Ji weaker than some cut-off | Ji| < Jcm are separated 
by a typical distance of order 

IcM ^ Jcm (4) 
The magnetization of the unfrustrated domain between two such weak bonds is of order 

mcM ^ y/lcM ~ JcaP (5) 



The flipping of such a domain thus involves a typical energy of order Jcm for the creation of two weak frustrated 
bonds, but allows to gain a 
optimal cut-off of order |^ 



bonds, but allows to gain a magnetic energy of order hmcM ^ ^^cm^ ■ "^^^ balance between the two terms yields an 



Jcm - h'/^'+^^ (6) 
so that the magnetization per spin Mg presents the following non-analytical behavior 

^ ^ ^ hT^ (7) 



CM 



The zero-temperature susceptibility 



^(T = 0) ^ ^ ^ /i-2/(2+M) (8) 



thus diverges at zero field /i — > 0. For instance, the Gaussian distribution p{J), which corresponds to the exponent 
/i = 1, leads to the behavior Alg ^ h^/^. The critical exponent (5+^ for the magnetization in external field was 
found to be exact via transfer matrix calculations by Gardner and Derrida Q, where the prefactor was moreover 
computed. The presence of small couplings does induce interesting new properties for the ground state with respect 
to the bimodal case. 

Chen and Ma have also analyzed the low temperature properties, in particular in the regime where the temperature 
T is much smaller than the typical energy Jcm of a domain. In this regime T <C Jcm ®j only a small fraction of 
two-level excitations will be excited : the density p{E ~ 0) of excitations near zero energy can be estimated to scale 
as 

= 0) ^ X ^ J^ ^ ;,2(,-i)/(2+,) (9) 
i-CM Jcm 

where 1/lcM represents the density of frustrated bonds in the ground state, and where 1/ Jcm represents the fraction 
of these frustrated bonds that will be involved in excitations of vanishing energy E 0. This excitations density is 
then expected to govern the leading term of the specific heat at low temperature 

C ~ T—p(E = 0) + O(T^) (10) 

T^O 6 

In this paper, our aim is to study in some details the statistical properties of the ground state configuration, and 
of the the low-energy excitations that govern the low-temperature properties. 



II. STATISTICAL PROPERTIES OF THE GROUND STATE CONFIGURATION 



In this section, we present the numerical results that we have obtained via the zero-temperature transfer-matrix 
formulation from which one can obtain the ground state configuration {at} in each given sample {Ji}. In the 
whole paper, we have used the following distribution for the couplings 

p(J) = ^iJjA'-i for -1<J<1 (11) 

yielding A = ^. The results given in this section have been obtained from averages over 10^ independent chains 
containing N ^ 10^ sites. 
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FIG. 1: Magnetization per spin in the ground state as a function of the external field varying between h = 5.10~* and 
h = 1.10~* in log-log plot, for disorder distributions IIH with exponents fi = 1,2,3. The corresponding slopes are in full 
agreement with the exact exponents 1/3, 1/2, 3/5 ( Eall2|l. 



A. Magnetization per spin 



As a first observable, we have computed the magnetization per spin which corresponds to a thermodynamic quantity 
which is exactly known from a transfer matrix calculation done by Gardner and Derrida who have obtained |5| 

^exact _ ^ ^) ( \ Y+2) _^ ^^^^ 



h^O \^(/X + 2) 



We have checked that both the scaling in h (see Figure ^ and the prefactor are in excellent agreement with the exact 
result l(T^. 

Now that we have identified the regime in h where the scaling ^ is satisfied, we may turn to more refined statistical 
properties for which, to the best of our knowledge, no exact expression is available. 



B. Probability distribution of frustrated links 



We have computed the normalized probability density of coupling |J| among frustrated links 

Pf(j) ^/('^) (13) 

where Nf{J) represents the number of frustrated links of strength J. 

According to the Chen-Ma argument, the frustrated links should have a typical strength of order Jcm ^ 
We have thus plotted on Figure [3 the probability distribution of frustrated links in terms the appropriate rescaled 
variable 

K=^^ = ^^ (14) 

JcM{h) /l2/(2+M) ^ ' 

for various h with the same initial distribution (|ll|l corresponding to /i = 1. 

To compare the distributions of frustrated links corresponding to different disorder distributions p{J) characterized 
by different exponents /i Hll(l . it is more convenient to consider the reduced variable 

r = |jr (15) 
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FIG. 2: Rescaled probability distributions P^^i{K = \ J\/ JcM{h)) of frustrated links in the ground state : the data for various 
fields, namely h — 1.10"'^ (line), h — 5.10~* (circles) and h = 1.10"'* (triangles), follow the same master curve. 




FIG. 3: Probability distributions P/[{R = r/rcAi{h) = {\J\/JcM{h))'^) of the frustrated links for various disorder distributions 
mi l corresponding to the exponents fJ. — 1 (bold line), /i — 0.5 (thin line) and ^ = 0.1 (dashed line). 



which is distributed with the flat distribution 

n''P"°"{r) - 61(0 < r < 1) 
for any /i in Taking into account the Chen-Ma scaling 

rcAiih) = J^.jih) = 



(16) 



(17) 



we have plotted on FigureOlthe probability distribution P^{R = r/rcM{h)) of the frustrated links for various disorder 
distributions corresponding to the exponents fi — I, fi = 0.5 and fj, = 0.1 : for /i = 1, this distribution is rather smooth, 
whereas it becomes steeper as /i decays. In particular, for fj, = 0.1, it becomes close to a simple theta function 9{R < 1). 
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FIG. 4: Rescaled probability distributions ^^(A = l/lcM{h)) of the length between two frustrated bonds of the ground state, 
for various initial disorder distribution characterized by the exponent ^ = 1 (full line) /j, = 1/2 (dashed line), ^ = 1/4 (circles) 

This shows that in the hmit /i ^ 0, the properties of the ground state become simpler, and we will discuss this point 
in more details in Appendix 1X1 

C. Probability distribution of the lengths of unfrustrated intervals 

According to the Chen-Ma argument, the length I between two frustrated links has for typical scale 1cm ~ 
/j-2m/(2+m)^ Indeed, we obtain that the appropriate rescaled variable for the length of unfrustrated intervals is 



A: 



I 



lcM{h) 



^/j2p/(2+M) 



(18) 



as h varies. The probability distribution ^^(A) of the scaling variable A is plotted on Figure 0] for various fi. 

For A ^ 0, in contrast with the bimodal case where 'P{\ = 1/Iim) presents an essential singularity 0|, we obtain 
here power-law behavior near the origin 



7'^(A) oc A"(^) 



(19) 



with an exponent a(/i) that grows with /i (see Figure 0}. For instance for /i = 1, the best fit yields the exponent 
a{n = 1) ~ 0.8. For large A, the decay is exponential 



^p(A) « e 



-7(m)A 



(20) 



For ^1= I, the best fit yields 7(^ = 1) ~ 0.25. 



D. Probability distribution of the magnetizations of unfrustrated intervals 

Similarly, we find that, in agreement with the Chen-Ma argument, the appropriate rescaled variable for the mag- 
netization of unfrustrated intervals is 



mcM{h) 



(21) 



as h varies. The probability distribution (5^(a;) of the scaling variable uj is plotted on Figure [S] for various /i. 
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FIG. 5: Rescaled probability distributions Qi_i{uJ — m/mcM{h)) of the magnetization between two frustrated bonds of the 
ground state, for various initial disorder distribution characterized by the exponent ^j. — 1 (fuU line) /i = 1/2 (dashed line), 
fj. = 1/4 (circles). 



Again, for w — > 0, we obtain power-law behaviors 



)^(w) cx w 



/3(m) 



with for instance /3(/i = 1) ~ 1.5, whereas the decay for large uj is exponential 
with S{fi = 1) ~ 0.5. 



(22) 



(23) 



III. LOW-TEMPERATURE PROPERTIES 



Low temperature transfer-matrix results for various observables 



We have first computed via transfer matrix various observables in the low-temperature regime T <^ JcM{h) 10 • For 
instance for h = 0.02 corresponding to JcM{h) = 0.07, we have checked, for temperatures T = 2.10~^, T = 3.10"^, 
T = 4.10^'^, T = 5.10^'^, that the leading term of the specific heat is linear in T 



C 



< El > - < En >^ 



oc T 



and that the leading term of the susceptibility is a constant 



X{T) 



< M% > - <Mn >^ 



TN 



■ oc cte 

T->0 



(24) 



(25) 



We have also found that the Edwards-Anderson order parameter deviates linearly in temperature from the zero- 
temperature value qEA{T = 0) = 1 



qEA = <<J^>^ ~ 1 - r(cte') 

T-+0 



(26) 



We have then studied the dependence in the external field h at fixed temperature. For instance, for T = 5.10 ^, 
we have studied the dependence in h = 10~^ to h = 10~^ (see Figure IHJ. The results for the specific heat is in good 
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FIG. 6: Low temperature behaviors of the specific heat C, of the Edwards- Anderson order parameter q and of the susceptibility 
X as a function of the external field in log-log plot for the special value fi = 1 : the exponents are respectively —0.04 for the 
specific heat, —0.66 for q and —0.6 for x (see text for more details). 



agreement with the Chen-Ma prediction (Eg 1101 and 1^ 

C 



/j2(m-1)/(2+m) 



as well as the susceptibility 



T T^O 



X{T) ~ /l-2/(2+A') 
T-*0 



(27) 



(28) 



We have also computed the Edwards- Anderson order parameter, and found the same exponent as for the suscepti- 
bility lEHl) 



1 - qsA 



^-2/(2+/.) 



T T^O 

which can be explained from the analysis in terms of the low-energy two-level excitations, as we now explain. 



(29) 



B. Interpretation in terms of low-energy two-level excitations 

We have already given the expression 11U|) of the specific heat in terms of the density p{E — 0) of excitations near 
zero energy. The similar expressions for the Edwards- Anderson order parameter and the susceptibility read Q 



1 ~ qEA 
T 



T~*0 



-\-oo 



dl lp{E = 0,0 



(30) 



/+CO 
dm ni^p{E = 0, m) 
-oo 



(31) 



where p{E — 0,1) represents the density of excitations of length I, and where p{E — 0, m) represents the density of 
excitations of magnetization m. Note that in the random field chain Q, the magnetization of a ferromagnetic domain 
is equal to its length, whereas here in the spin-glass chain, it is not the case, since the magnetization of an unfrustrated 
domain scales as This is why the scaling in h are the same here for these two observables. 
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FIG. 7: Nature of the low-energy excitations. The ground state is made of long unfrustrated intervals separated by frustrated 
bonds called ..io, j'o, ko, lo---- Low energy excitations may be of several types (displacement, annihilation or creation of frustrated 
bonds). 



We have numerically computed the probability density of excitations, as a function of their size I and of their type. 
Indeed, we obtain by an exhaustive numerical enumeration that the total density of low-energy excitations is exactly 
the sum of three contributions 



ptot{E = 0, - pli%{E - 0, + p'aLiE = 0, + £rL{E = 0, (32) 

(1) The excitations of type 1, of density pl^^gpiE = 0,1), involve a single frustrated hnk of the ground state, 
of coupling Ja, that can be displaced to another position of coupling Jb (inside the intervals defined by the two 
frustrated neighbors of Ja) with almost no energy cost. The energy difference 

AE'-^^ ~2\Ja\ + 2\Jb\ + 2hmab < T (33) 

involves the defrustration of Ja, the frustration of Jf, and the magnetization niab between these two links. 

(2) 

(2) The excitations of type 2 involves a pair of neighbor frustrated bonds that can be annihilated {Pannii^ — 0; 0) 

f 2") 

or that can be created {pcrea{E — 0, /)) with almost no energy cost 

A^l, = -2| Jil - 2| J2I + 2hm,2 < T (34) 
AE^la = 2| Jil + 2| J2I + 2/imi2 < T (35) 

These two type of excitations are symmetric and are thus expected to correspond to the same distribution 

PaL{E^O,l)^p^la{E^O,l) (36) 

which we have checked in our numerical results. The densities of these excitations in terms of their rescaled length 
A — l/lcM{h) are given on Figure |H1 for ^ = 1. 

We have moreover checked the relations (|1UI3UI31|I between, on the one hand, the specific heat, the Edwards- 
Anderson order parameter, and the susceptibility obtained from low temperature transfer-matrix calculations, and 
on the other hand the total number of excitations, their averaged length, and their averaged square magnetization. 
The agreement shows that the excitations described above are the only ones that play a role in the low temperature 
behavior of these observables. The present analysis in terms of the statistics of low-energy excitations thus gives a 
microscopic interpretation of the low temperature equilibrium properties. 
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FIG. 8: Densities of the three types of low energy excitations present in equation l|32|l . in terms of their rescaled length 
A — l/lcKiih) in the case /i = 1. 



IV. 



RENORMALIZATION PROCEDURES IN EACH SAMPLE 



As already mentioned in the introduction, in the bimodal case ±J, there exists a real-space RG procedure that 
allows to compute exactly the statistical properties of the ground state (such as the domain length distribution) 
as well as the statistics of low-energy excitations [3] . A natural question is thus : is there a generalized RG procedure 
that would be valid for the spin-glass beyond the bimodal case Ji = ± J? Before trying to answer this question for 
the p( J) ~ ^1 J|'^~^ case, let us first briefly recall the principle of the RG procedure for the bimodal case. 



A. Bimodal distribution 

The RG procedure defined in ref ^ for the bimodal case, consists in an optimization from small scales towards 
large scales : one starts for instance from the completely magnetized state ct^ = -1-1, that contains many frustrated 
bonds, and one flips iteratively the unfrustrated domain presenting the smallest magnetization, as long as the energy 
is lowered, i. e. as long as the balance between the energy gained from the suppression of the two boundary frustrated 
bonds is bigger than the energy loss {2hm) from the negative orientation with respect to the external field 

AEf^'P = -4 J + 2h\m\ < (37) 

So the RG procedure has to be stopped when all unfrustrated domains have magnetizations m > 2J/h : the state 
obtained then corresponds to the ground state. What makes the renormalization tractable in this case is that due 
to the constant cost (2 J) of any domain wall, the renormalization concerns a one-dimensional potential, namely the 
magnetization as a function of the running point. 



B. General distribution : hierarchical RG based on the energy 

In this case, the problem cannot be reformulated as the renormalization of a one-dimensional potential, since in 
addition to the magnetization, one has to take into account that the couplings Ji vary along the chain. Moreover, we 
have seen with the Chen-Ma argument that the frustrated links are concentrated on small couplings. 

As a first step, we can thus formulate the following renormalization that optimizes from the biggest scales towards 
smaller scales : 

(i) one starts from the state with no frustrated links that presents a positive magnetization (of order ^/N for a 
chain of size TV) : it corresponds to one of the two mirror zero-field ground states. 
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(ii) at each step, we flip the interval (cTi+i, ...aj) that allows the maximal decrease of energy E = min AE{i, j) < 

i<j 

where 

i 

AEii,j)^2f{J,)\M + 2fiJ,)\J,\ + 2h (38) 

k=i+l 

where f{Ji) = 1 if Ji becomes frustrated during the flip, /(J^) = — 1 if becomes unfrustrated during the flip. If we 
find no interval to flip, we have obtained the ground state, since by definition the ground state is stable upon the flip 
of any interval. 

We have indeed checked that this procedure allows to obtain the exact ground state computed independently via the 
zero temperature transfer matrix. We have moreover obtained that actually, the links that become frustrated during 
the procedure never get 'unfrustrated' later in this procedure. And in fact the following hierarchical RG procedure, 
where at generation n, the chain is cut into a certain number of sub-chains, gives the exact ground state : 

(i) same initial state as before 

(ii ) once the first interval to be flipped is found, we can find the next intervals to be flipped independently 

within the three sub-chains (l,ii), and {ji,N). And we iterate until there is no interval to flip anymore. 

The fact that the first links that become frustrated indeed belongs to the final ground state, can be justified 

via a 'reductio ad absurdum' ,7], and then it is valid for all stages of the procedure. 

This hierarchical procedure thus defines an energy driven fragmentation process of the chain, whose statistical 
properties can be studied and compared to other fragmentation models In particular, we have studied the number 
nspHt of splitting and the number Ugene of generations as a function of the size N of the chain, for N — 500 to 

= 8000. For the disorder distribution (|11|) with fJ, = 1 and external field h = 0.02, that corresponds to the length 
Chen-Ma scale Icnih) ~ 13.5.., we obtain, as expected, that the number of splittings grows hnearly in N 

A^ N^oo 

and that the number of generations grows logarithmically in A^ 



lO"'' (39) 



In A^ AT^oo 



2.1 (40) 



This RG analysis reveals a hierarchical structure among the frustrated links of the ground state. This hierarchy 
has both a spatial meaning, but also an energy meaning. Indeed, since an interval created inside another interval 
has, by definition of the RG procedure, a smaller energy, it is clear that the low-energy excitations of the type 
'annihilation' introduced in 1)32(1 can only concern a pair of frustrated bonds that have been created together and that 
have no descendent in the hierarchy. Similarly, the low-energy excitations of the type 'creation' introduced in l|32|) 
correspond to a pair of frustrated bonds that would have been created next if the procedure had been applied a bit 
beyond AE < 0. Finally, the only remaining excitations are the 'displacements' H32|l . that actually also preserve the 
hierarchical structure, since a link can move only between its two frustrated neighbors. 



C. RG based on the weakest link at each step 



Since the frustrated links concentrate on the links that are weak, i.e. of order JcM{h), it is tempting to try to 
define a RG procedure based on the weakest link at each step. In Appendix^ we show that a simple RG procedure 
based on this idea becomes exact in the limit /i 0, that corresponds to an infinitely broad distribution 1111) near 
0. 



V. CONCLUSION 



For the spin-glass chain in an external field h, we have first studied via zero-temperature transfer-matrix calculations 
the statistical properties of the ground state configuration. We have then studied the nature and the statistics of the 
low-energy two-level excitations, via a direct enumeration, and we have analyzed their contribution to the specific 
heat, the Edwards-Anderson order parameter and the susceptibility in the low-temperature regime. Finally, we have 
shown that an extended RG procedure, based on the iterative flipping of the best energetic interval, could be used 
to obtain the exact ground state in external fleld. This RG procedure reveals a hierarchical structure among the 
frustrated links present in the ground state. 
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The possible relation of this hierarchical picture with higher dimensional disordered models is clearly of interest. In 
the 2d random field Ising model, a spatial hierarchical picture has been identified long ago , based on the existence 
of the Imry-Ma domain length scale. We tend to think that this hierarchy is energetic in nature. More precisely, we 
believe that a RG procedure that would start from the ferromagnetic pure state, and flip iteratively, at each step, 
the most advantageous domain- regardless of its size- will ultimately converge towards the ground state. After the 
initial flipping of the most advantageous domain (which, through the Imry-Ma argument, is also the largest), one has 
to search separately for the next advantageous domains inside and outside the initial one. These flipped domains will 
then display a disjoint or hierarchically nested structure. Another related problem, where an iterative optimization 
procedure starts from the largest scale and hierarchically proceeds towards smaller scales, has been studied by Binder 
[lO| for interfaces. 

For the 2d or 3d spin glass case at zero magnetic field, there is no equivalent of the Imry-Ma or Chen-Ma length 
scale. The existence of a hierarchical organization in ground state or low temperature properties has nevertheless been 
found along various lines: rigidity properties jTll. distance between spin configurations subject to the same thermal 
noise 0, 0| , calculations on small systems |l4j. extensive data clustering analysis |0| . This hierarchical organization 
pertains to spin clusters, and is not a priori linked to a hierarchy in their flipping energies |lfi| . 
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APPENDIX A: SIMPLE SOLUBLE RG PROCEDURE IN THE LIMIT ^ ^ 

In the limit /i 0, the disorder distribution p( J) Hll() becomes infinitely broad distribution near J — s- 0. It is thus 
tempting to define a RG procedure based on the weakest link at each step. In the following, we consider the simplified 
RG procedure : 

(i) same initial state as before : one starts from the state with no frustrated links that presents a positive magne- 
tization 

(ii) First iteration : we choose the smallest coupling in absolute value F = | JminI- The chain is thus decomposed 
into two sub-chains. We consider the magnetizations of the two sub-chains mi -|- m2 = M > : if one of the two 
magnetizations (mi, 7712) is negative, for instance, mi < 0, we will flip the sub-chain 1 if it lowers the energy, i.e. if 
the balance between the cost 2F of introducing a frustrated bond between the two sub-chains is less than the energy 
gained by the orientation of the sub-chain 1 along the external field 

AE^^'P = -2/i|mi| + 2r < (Al) 

otherwise, if AE^^'^p > 0, we do not flip the sub-chain 1. After this, the two sub-chains will evolve as two independent 
sub-chains with free boundary conditions, so we iterate the procedure. 

This very simple RG procedure is of course not exact, since at each step, we neglect the weak bond at the other 
boundary of the interval that has been previously decimated | J| < F. Indeed, at each step, we consider that the cost 
of the flipping of an interval is exactly 2F, whereas it should be 2(F ± e), where e < F is the absolute value of the 
coupling at the other boundary, that has been previously decimated, and where the sign (±) depends on the state of 
this coupling e, frustrated or not, in the renormalized chain at the RG scale F. However, we will show below that it 
becomes exact in the limit /i — > 0, where the distribution becomes infinitely broad distribution near J — s- 0. 



1. Statistical properties of the intervals between weak bonds at RG scale P 

In this section, we study the RG procedure defined above in the thermodynamic limit of an infinite chain. At the 
renormalization scale F, the chain is split into independent unfrustrated intervals, separated by weak bonds that can 
be either frustrated or not. We now derive some statistical properties of these intervals between weak bonds. 

The distribution Pr{l) of the length (j — i) between two weak bonds Ji_j < F is simply exponential, since it 
corresponds to the probability that / independent couplings have \Jk\ > F, 



(A2) 
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with the characteristic length 



Ir 



hi 



r-o 2Arf' 



(A3) 



The magnetization of an unfrustrated interval of length / is simply the sum m = ± X)i=i sgn{Ji). Since I is large, 
the distribution of m given I is a Gaussian in the ground state at ft. = 0, and thus we obtain after averaging with 
respect to the length I with (jA2|) the following a priori distribution 



jjaprtort 



(m) 



' „ 1 _ie-Tr 1 

dl — e 'r = ( 

V2Td 2mr 



with the characteristic magnetization 



TOr 



1 



a) 



(A4) 



(A5) 



2 r^o 2r^ 

Now from this a priori distribution that describes the magnetizations of these domains in the h = ground state, 
we wish to compute the distribution of domain magnetizations obtained via the renormalization procedure, where we 
have tried to flip intervals in a iterative way. 

A domain existing at scale F was created at some previous scale T' representing the biggest of the couplings at the 
two boundaries. Since the distribution of the already decimated coupling reads 



Oi\J\<T)p{\J\) 



jdje{\j\<r)p{\j\) 

the distribution of the creation scale F' is simply 



9i\J\<Ty 



^creat^on^p/-) ^ 2pf'''°"(F') / d| j' jp^"'"" (| j' | ) ^ 6I(F' < F) 

^|j'i<r' 

r' 



(A6) 



(A7) 



The stability condition m > — at the scale F' of its creation immediately yields the simple properties for the 
probability distribution Pr{m) of the domain magnetization at scale F 



m > ^ ) Pr{m) = 2P^P"°"{m) 



-- > m ] Prim) = 



(A8) 
(A9) 



For the values \m\ < T/h, the probability distribution [Pr(™)] 
stable at the creation scale F' reads 



stabi 



F 



e[l.>m>Oj [Pr (to)] = Pr''"''"im) J dV 



induced by the only condition to have been 
,2^(F')2^-i 



papriori 



(to) 



1 



F' 



(AlO) 
(All) 



> m > -- ) [Pr(TO 



p2M 



papriori 



(to) 



1 



Imlft 



2/^ 



F' 

> TO > - — 



(A12) 
(A13) 



In the limit /i — > 0, we have thus the simplification that negative m become negligible, because two bonds weaker 
than F are typically much weaker than F, as a consequence of the broadness of distribution. So at leading order in /z, 
we have the simple property 



[Pr(TO)] 



stable 



(m > 0)2P^P"°"{' 
1 

(to > 0) e 

TOr 



(A14) 
(A15) 
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2. Probability measure for the fragmentation process at scale F 

The probability measure to find a bond of strength F inside an interval (L, M) existing at scale F" that becomes 
fragmented into two sub- intervals (Zi,mi) and (Z2,'Ti2) reads 

Mt{L, M; li,mi;l2, m2)d^dLdl\dl2dm\dm2 (A-16) 

= drpr{r) — e~'^dkdl26{L - {k + l2))dM20{M > 0)dmi-^dm2-^5{M - (mi + ma)) (A17) 
It v27rZi \/1-k12 

Indeed, we have the following properties, for the integration over some variables. After the integration over (mi, ma), 
the distribution of M is a Gaussian as it should 

J dmi J dm2Afr{L,M;li,mi;l2,m2) = (A18) 

m2 

pr(r)-e-*5(L - ih + l2))2e{M > 0)^^ (A19) 

After the integration over all magnetizations {mi, 1712. M), the distribution for (^1,^2) is uniform except for the con- 
straint I1 + I2 = L 

j dM j dmi j dm2Afr{L,M;li,mi;l2,m2) = (A20) 

pr{T)^e-^5{L-{h+l2)) (A21) 
tr 

After the integration over (^i, Z2), the probability to find a bond F in an interval of length L 

j dh j dh j '^^ j '^"^i j dm2Mr{L,M;h,mi;l2,m2) = pr{T)^e~^L (A22) 

is proportional to LPy{L) since they are L possible positions. 

In the following, when computing observables concerning the flips at scale F, it will be more convenient to integrate 
first over the lengths that play no direct role in the flip condition, to keep the magnetizations that enters the flip 
condition 

M'(M;mi,m2)= j '^^ J '^^i J c^^2A/'r(L, M; Zi, mi; ^2, TO2) (A23) 

= Pr{r) — 20{M>Q)6{M-{mi+m2)) dhe''^ d^e'^^-^ {A24) 
Ir Jo \/27rZi V^nh 

= priT)e{M > 0)S{M - (mi + m2))e""'^e""'^ (A25) 

3. Flipping probability at scale T 

The number of bonds of strength F that becomes frustrated at scale F is proportional to 



Af/™'(F) = jdMjdmij dm2Afr{M;mi,m2) 



{mi< ~] + 6 (m2 < 



h \ h 



(A26) 



r f i^ii iM-mii / r\ 

pr(F) / dM0{M > 0) / dmie'^e ^^26* ( "^i < "-^ j (A27) 

2/9r(F) / dMe ^ / dml^e ^ (A28) 

pv{T)mle~'^ (A29) 
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that should be compared with the total number of bonds of strength F proportional to the normalization 



iV^°*(r) = / dAf / dmi I dm2AfriM;mi,m2) 



r+oo 

Pr(r) / dM 



dmie '"r e 



"r + / dmie "re 
/o 



+ 00 



dniie 



M 



= 2pr(r)m2 

The flipping probability Fr (F) of a bond T at scale F is given by the ratio of the two 



where the rescaled variable 



2 /A 



1/2 



(A30) 

(A31) 
(A32) 

(ASS) 
(AS4) 



represents the ratio between the quantity T/h appearing in the flip condition and the typical scale mr of the magne- 
tization of a domain existing at scale F. At the beginning of the procedure x <^ 1, there is a finite probability of flip, 
of order 1/2, whereas for x ^ 1, the probability of flip becomes exponentially small. 



4. Magnetization per spin at the end of the procedure 



To compute the magnetization per spin TOspm(^)j we have to integrate over all the flips done at various scales, and 
to keep track of the associated magnetization gain 



P Si J 



(Am) 



r,r+dr 



Ir 



(A35) 



where (Am)r,r+(ir is the mean magnetization gain associated to a domain flip at scale F, which can be expressed in 
terms of the measure (IA25I) 



(Am)r,r+dr = / dM / dmi / dm27Vr(Af; toi, 7712) 



pr(r) / dMOiM > 0) / dmie 



nr g 



mi < j 2|mi| + e (m2 < ) 2|m2| 



F 



nil < --^ ) 4|mi| 



r+°° M 

= 4pr(F) / dMe ~ 
= pr(r)m^(l + 2a;)e-2^ 



anil m^e 



(A36) 

(A37) 

(ASS) 
(A39) 



where x is the scaling variable defined in ljA34p . 

Finally, using pr(r) = 2AT^^^ and the new variable x as integration variable instead of F, we obtain the magneti- 
zation per spin (jA35p as 



r+oo 

{h) = (1/2) / dFpr(r)mr [1 + 2x] e 
Jo 



-2x 



(A40) 



and the result 



mspin{h) = /if +2 ( — 1 CsimpleifJ-) 



(A41) 
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The exponents in h and A agree with the exact results of whereas the prefactor reads 



^simple 



2(/i + 2) 



Jo 2 + V 



1 , 1 - lEuler - In 2 2 

2+/^ 4 



(7iJ«ier+ln2)2 12 _ ^2 
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instead of the exact prefactor obtained via transfer matrix computations 5] that reads in our notations 
^ ^ / , o\ — , 1^ 1 l-7E„,er-ln2 2 (Tfiufer + In 2)^ 3 

Cexactif^) = (M + 2) + 1) / ^ X = 2 + ^ 4 + 16 

^ ' 11+2 



(A42) 
(A43) 



(A44) 



so the discrepancy with the exact prefactor (|A44p only appears at order ix^. 
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